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ABSTRACT 

There is a growing body of evidence for the presence of multiple stellar populations 
in some globular clusters, including NGC 1851. For most of these peculiar globular clus- 
ters, however, the evidence for the multiple red giant-branches (RGBs) having different 
heavy elemental abundances as observed in co Centauri is hitherto lacking, although 
spreads in some lighter elements are reported. It is therefore not clear whether they 
also share the suggested dwarf galaxy origin of co Cen or not. Here we show from the 
CTIO 4m UVI photometry of the globular cluster NGC 1851 that its RGB is clearly 
split into two in the U — I color. The two distinct RGB populations are also clearly 
separated in the abundance of heavy elements as traced by Calcium, suggesting that 
the type II supernovae enrichment is also responsible, in addition to the pollutions of 
lighter elements by intermediate mass asymptotic giant branch stars or fast-rotating 
massive stars. The RGB split, however, is not shown in the V — I color, as indicated by 
previous observations. Our stellar population models show that this and the presence 
of bimodal horizontal-branch distribution in NGC 1851 can be naturally reproduced if 
the metal-rich second generation stars are also enhanced in helium. 

Subject headings: globular clusters: individual (NGC 1851) — globular clusters: for- 
mation — stars: abundances — stars: evolution — stars: horizontal-branch 



1. INTRODUCTION 

The first photometric evidence for the presence of multiple stellar populations in globular 
clusters came from the discoveries of the discrete distributions of stars on the r ed giant-branch 



(RGB) and main -sequence (MS) of the most luminous globular cluster u Cen (|Lee et al 



1999 



Bedin el al.1 12004). The fact that the distribution of RGB stars is not just showing a spread but is 
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discrete is a compelling evidence for the heavy elements enrichment and the formation of successive 
metal-enhanced generations of stars in proto-w Cen. The similarity of this feature to that of the 
M54, suggested as a nuclear star cluster of the Sagittarius dwarf galaxy (jLavden fc Saraiedini 
19971 ). together with other peculiarities of ui Cen, has led the community to conclude that to Cen 
was once part of a more massive dwarf galaxy that merged w i th the Milky way, as the Sagittarius 



dwarf galaxy is in the process of doing now (jLee et al 



1999 



Bekki Freeman! 120031 ) . In recent 



years, more and more evidence is reported for the presence of dou ble or multipl e stellar population s 



in other glo bular clusters, such as NGC 280 8 (jPiotto et a" 



NGC 6388 (jMoretti et alJ 120091 ) . and M22 (jMarino et al 



iOOa- iDa ( os(a et al 



2007 ) NGC 1851 dMilone et alJl200«h 



20091 ). For most 



of these peculiar globular clusters, however, the evidence for the discrete distribution of heavy 
elements as observed in the R GB of u Cen is hitherto lacking, a lthough spreads in some lighter 



2005 



Piotto et al 



2007: 



of M22 (Marino et al 



elements (ICarretta et al.ll2008l. and refere nces therein) and helium (jLee et al 



Yoon et al 



2009 



sp 

l.: 



2005 



D'Antona et al 



2008) a re rep orted. Therefore, the case of u Cen, and now that 



Da Costa et al.ll2009l ). are still viewed as exceptional, and the presence 



of chemical inhomogeneity in other globular clusters is largely considered as due to the pollution 
mechanisms expected in normal globular clusters, such as the winds from the intermediate mass 
asymptotic giant bra nch (AGB) stars or fast-rotating massive stars (| Ventura &: D'Antonal 120081 ; 
Decressin et al J 120071 ). 



The purpose of this Letter is to report that one of these globular clusters, NGC 1851, is showing 
a clear split in the RGB from our C TIO 4m UVI photometry. This is compared with the Ca-by 
photometry (IJ.-W. Lee et al.ll2009al ) to confirm that the two distinct RGB populations are different 
in the abundance of heavy elements. Stellar population models are then constructed to show that, 
when the metal-rich second generation stars are also enhanced in helium abundance, the split of 
the RGB discovered in U — I color would not be detected in usual optical colors, such as V — I 
used in the HST/ACS survey. 



2. OBSERVATIONS AND COLOR-MAGNITUDE DIAGRAMS 



Our observations were performed using the CTIO 4m Blanco telescope during 2007 November 
13—16 and 2008 October 27—31. The telescope was equipped with the Mosaic II CCD Imager 
consisting of eight 2k x 4k SITe CCDs providing a plate scale of 0.27 arcsec pixel -1 and a field of 
view of 36 x 36 arcmin 2 . All of our science frames were obtained under photometric conditions. The 
total exposure times for UVI were 3990, 676, and 573s, respectively, split into short, intermediate, 
and long exposures in each band. NGC 1851 was placed on chip 6, a pproximately 4.5 arc min South 



and 6 arcmin East from the CCD center. Several standard fields ( Landolt 



1992 



20071 ) were also 



observed during our observing runs. The IRAF0 MSCRED package was used for preprocessing 



X IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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including bias correction an d flat fielding. PSF photometry was then carried out u sing DAOPH OT 
II/ALLSTAR(|Stetso3[l983), and D AO GROW was used for aperture corrections (jstetsonlll990r ). 



Figure 1 shows color-magnitude diagrams (CMDs) of NGC 1851 in (U, U — I) and (V, V — I) 
planes. Stars within 3.6 arcsec and outside of 4.5 arcmin from the cluster center are excluded from 
the CMDs to reduce blending effects and the field star contaminati ons, respectively. To examine the 
CMD features more carefully, we adopted the "separation index" (jStetson et al.ll2003l ) for selecting 
stars that are relatively less affected by adjacent starlights. All the stars in our CMDs lie well 
within chip 6, and therefore our CMDs are not subject to any uncertainty stemming from the 
possible chip to chip variations of the m osaic CCDs. O pen circles denote RR Lyrae variables in our 
program field among those identified by IWalkerl (|1998l ). plotted at random phase of pulsation. The 
most remarkable feature of Figure 1 is the presence of two distinct RGBs in the (U, U — I) CMD. 
The discrete distribution is clear from the sub giant-branch (SGB) to the tip of RGB where the 
mean separation on the RGB is ~0.27 mag in U — I. When measured at given I magnitude, this 
value is red uced to ^0.20 mag. For the bright RGB stars, some early hints for t his feature were 



noted b y ICalamida et al.1 (|2007l ) from the Stromgren (ml, u — y) CMD , and by 



J.-W. Lee et al 



( 2009a ) from the Ca-by photometry (e.g. see their Figure 1). Note that Milone et al. ( 20081 ) only 



discovered a split in the SGB, and the split of RGB we discovered in this paper was not detected in 
their HST/ACS photometry employing F606W and F814W passbands. Similarly, the RGB split 
is not apparent in our (V, V — I) CMD (Figure lb). This is most likely because the [/-band is more 
sensitive to metal abundance variation than other passbands, as more metal atomic and molecular 
lines are located in the J7-band (see section 3). 

Given the small foreground reddening value of E(B— V) = 0.02 ( Harrisfl996 ) toward NGC 1851, 
it is very unlikely that the differential reddening has caused the double RGBs. The color difference 
between the two RGBs in U — I color, at given / mag of the horizontal-branch (HB) level, is 
~0.20 mag, which is about three times larger than the maximum color difference expected in the 
extreme situation where one group of stars are all reddened by E(5 — V) = 0.02, while the other 
group has E(£> — V) = 0.00. Also, if the observed color difference in U — I color is due to the 
differential reddening, it would result in the V — I color difference of 0.10 mag, which would have 
been detected along the RGB in our (V, V — I) CMD. Similar spatial distributions of stars on the 
bluer and redder RGBs also indicate that the differential reddening, if any, is not likely the cause 
of the double RGBs. 



3. DISCUSSION 

We have shown that the RGB of NGC 1851 is split into two distinct subpopulations. The 
fact that the distribution of RGB stars does not just show a spread but is discrete can naturally 
eliminate the possibilities such as (1) star formation from not well mixed inhomogenous interstellar 
matter, (2) differential reddening (see also section 2), and (3) photometric errors, as all of these 
would produce a spread in color rather than a discrete distribution. Consequently, the most 
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plausible interpretation for the double RGBs is that the NGC 1851 underwent metal enrichment 
in its early stage of evolution and successively formed metal-enhanced second generation stars. 
Spectroscopic observations show star-to-star abundance vari ations of the li ghter elements (elements 
lighter than Si such as N, O, Na, & Al) in NGC 1851 (|Hesserl Il982l : lYong fc Grundahll hood : 
Yong et al.ll2009r ) . This is generally interpret ed as a result of pollution from winds of intermediate - 
mass AGB stars ((Ventura &; D'Antonall2008l ) or f ast-rotating massive stars (IDecressin et al.l l2007'l. 
Based on the elemental abundances observed by lYong &; Grundahll (|2008l ) and lYong et al.1 ((2009) 
for eight bright RGB stars, we can estimate t he differences in t he ligh ter elements between the two 
RGB populations. Analysis of these data by |j.-W. Lee et al.l (|2009al . see their Figure 3) indicates 
that, while N, Na, and Al are all enhanced in redder RGB (A[N/Fe] w 0.47, A[Na/Fe] « 0.53, 
A[A1/Fe] ~ 0.20), O is depleted (A[0/Fe] ~ —0.45), and no significant variation is shown in 
[C/Fe]. 

In order to investigate the effect of these elemental variations in the (U — I) color, in Figure 
2, we have computed the differe nces in fluxes between t wo synthetic spectra using the spectral- 
synthesis program SPECTRUM (IGray &; Corballyl Il994l ). one without and the other with these 
enhancements (and depletion for O) in lighter elements (panels a & b). Also shown in Figure 2 are 
for the two additional cases where (1) CNO and Na are enhanced by 0.3 dex while heavier elements 
are fixed (panels c & d), and (2) elements heavier than Al are enhanced by 0.3 dex while lighter 
elements are fixed (panels e & f). These simulations demonstrate that small variations either in 
lighter or heavier elements could cause significant change in (U — I) color. The observed variations 
in the lighter elements alone, however, would cause relatively small line blanketing equivalent to 
A(U — I) ~ 0.079 in terms of color difference. Despite uncertainties, taken at face value, this is 
only ~41% of the observed color difference [A(U — I) = 0.195 ± 0.011] at given I magnitude and 
suggests that other effects may also be responsib le for the RGB split in (U — I) color. Indeed, 
recent Ca-by photometry (IJ.-W. Lee et al.ll2009al ) shows that besides lighter elements variations, 
RGB stars of NGC 1851 also show bimodal distribution in Ca, which can only be supplied by Type 
II supernovae (SNe II; iTimmes et al.l Il995j). In Figure 3, we sho w in our (U, U — I) CMD the 



"Ca- normal" and "Ca-strong" stars from IJ.-W. Lee et al.l (|2009al ). The Ca-strong stars lie well 



on the redder RGB seq uence, whereas the Ca-normal stars are on the bluer RGB. According to 
J.-W. Lee et al.l (|2009bl ). the difference in Ca abundance is estimated to be A[Ca/H] « 0.15 dex, 
and other heavy elements, albeit small, are similarly enhanced in redder RGB population. We 
conclude therefore that the redder RGB population is richer than the bluer RGB population not 
only in lighter elements (N, Na, & Al) but also in heavy elements (such as Ca, Si, Ti, & Fe). 

In order to better understand the origin of the RGB color difference in (U — I), and to place 
stronger constraints on the chemical combinations of the two subpopulations, we have constructed 
stellar population models based on the updated version of the Yonsei-Yale (Y 2 ) isochrones (Yi et 
al. , in prepara t ion) a nd HB evolutionary tracks (Han et al., in preparation). Readers are referred 
to lYoon et al.l (|2008l ) and references therein for the details of our model construction. Figure 4 
presents our synthetic CMDs for NGC 1851 in (U, U — I) and (V, V — I) planes. Our models 
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are constructed under two different assumptions regarding the chemical enrichment in NGC 1851. 
First, we assumed that the second generation population is more enh anced in met a llicity , but 
not in helium (Fig ure 4a; hereafter AZ-only model). Note that, while ICassisi et al.l (|2008l ) and 
Salaris et al.l (|2008l ) suggested the difference in only CNO abundance, here we are assuming the 
difference in overall metallicity of 0.15 dex as discussed above. Second, in Figure 4b, we then 
assumed that both metal and helium abundances are enhanced (hereafter AZ+AY model). The U 
flux is very sensitively affected by CN and NH bands, and therefore by N abundance. In order to 
reflect this effect in the (U, U — I) CMD of the second generation population (redder RGB), we are 
also including the effects of the additional line absorptions from the enhancements in N and other 
lighter elements as discussed above, again by using SPECTRUM. More rigorous modeling should 
include the effects of the se lighter elements e nhancements in the construction of stellar evolutionary 
tracks, but we note that iDotter et al.l (|2007l ) found the enhancement in N has only negligible effect 
in the HR diagram morphology at globular cluster ages. 

Comparison of the number ratio between the two subpopulations suggests that bluer RGB, 
brighter SGB, and redder HB are associated with one subpopulation, which we refer to as "Pop-1", 
while the redder RGB, fainter SGB, bluer HB are associated with the other subpopulation, which 
we refer to as "Pop-2". Pop-1 comprises about 75% of the total population, while Pop-2 takes 
about 25% of the whole populatiorJE The AZ-only model in Figure 4a (AZ = 0.0004 and Aage = 
0.1 Gyr) matches well with the observed CMD from the MS through the RGB in (U, U - 1) CMD. 
Yet, the model fails to reproduce the blue HB and the narrow RGB in (V, V — I) CMD (Figure 
4a, inset). This is because metal enhancement i n Pop-2 moves HB to the opposite direction, i.e., 
makes HB morphology redder (see for example iLee et al.lll994l ). At the same time, the increase 
in metallicity, albeit small, w ill create a gap in (V — I) color between the two RGBs which is not 
shown in the observed CMD (jMilone et al.ll2008l ). Therefore, the AZ-only model is in conflict with 
the observed CMDs of NGC 1851. Figure 4b shows that the AY+AZ model (AZ = 0.0004, Aage 
= O.lGyr, and AY = 0.05) is in good agreement with the ob servation from the MS to the HB. 
Note that, in our modeling, we are just employing the standard iReimersI (|19771 ) mass-loss law, and 
the same mass-loss parameter r) for both normal helium and enhanced helium populations. Input 
parameters adopted in our AY+AZ model are listed i n Table 1. The enhanced helium in Pop- 
2 easily overcomes the effect of metallicity on HB (see ILee et al.ll2005l ). shifting HB morphology 
toward blue and reproducing the observed HB bimodality. At the same time, the increase in helium 
abundance in Pop-2 moves RGB slightly bluer and this effect cancels out the metallicity effect in 
(V, V — I) CMD, making apparently single and n arrow RGB. Also, the model (Figure 4b, inset) 
is in better agreement with the observed SGB split ljMilone et al.ll2008l ). The [/-band is much more 
sensitive to metal line blanketing, and therefore the increased helium has relatively small effect on 
U — I color of RGB. One caveat to this helium enhanced scenario is that the blue HB in our model 



2 We have compared the ratio of the two subpopulations with samples selected in different manners (e.g. entire 
sample and a sample selected with various separation indices) and for all cases, the ratio co mes out to be s imilar . 
This is also more or less consistent with the population ratio based on the SGB stars reported in iMilone et all ([2009). 
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appears to be slightly brighter (~0.05 mag) than the observation in [/-band, although this could 
be due to the uncertainty in bolometric correction of the [/-band. 

The apparent helium enhancement in the second generation subpopulation in NGC 1851 is 
reminiscent of the cases of other globular clu sters with multiple popula t ions, including uj Ce n 



jNorrislhooilLee et al.lhooj: 



Piotto et al 



NGC 6388 and NGC 6441 jCaloi k D'AntonalhoO?! : lYoon et alJhood ). Although the origin of this 



2005]), NGC 2808 dLee et al 11200.4 iD'Antona et al.ll2005h 



helium enhancem ent is curren t ly not fully unders tood, general consensus is that it can be supplied 



either by SNe II (jNorris 



2004 



Piotto et al.ll2005l) or by the intermediat e mass AGB stars or fast 



rotating massive stars ([Ventura &: D'Antonall2009l ; lDecressin et al.ll2007l ). In the case of NGC 1851, 
given the enhancements both in the lighter (such as N) and heavy (such as Ca) elements in the 
second generation population (Pop-2), all of the mechanisms seem to be responsible for the helium 
enhancement. A possible scenario is that, soon after the formation of the first generation stars 
(Pop-1), numerous SNe II explosions enriched both metal and helium of the leftover gas in the 
proto-NGC 1851. Winds and ejecta from the intermediate mass AGB stars may have added more 
helium and simultaneously enhance lighter elements. The second generation stars (Pop-2) would 
have then formed from the gas enr iched in overall metallic ity, helium, and lighter elements. We note 
that the present mass (~1O 6 M ; iPrvor Sz Mevlanl Il993) of NGC 1851 is too small to retain the 



ejecta from n umerous SN explosions ( Baumgard t et al. 



2008 



Dopita Smith! Il986l ) . Therefore , 



as in uj Cen (|Lee et al.lll999l : iBekki Freeman! 120031 ) . we suggest that NGC 1851 is most likely 
the remaining core of more massive prim eval dwarf galax ies that merged and disrupted to form 
the proto-Galaxy, as recently proposed by lLee et al.l (|2007l ) for the origin of globular clusters with 
the extended HB. Our result presented here is calling a new survey of globular clusters employing 
UV and/or Calcium filters in order to detect small difference in metallicity expected in the yet to 
be identified globular clusters with multiple populations. High resolution spectroscopy of stars in 
the two distinct groups are also needed to confirm the small difference in the abundance of heavy 
elements. 
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Fig. 1. — Color-magnitude diagrams for NGC 1851. In both panels, only the stars within an 
annulus between 3.6 arcsec and 4.5 arcmin from the cluster center and with sep > 1.0 have been 
plotted. Note the discrete double RGBs in the U vs. U — I CMD. Open circles denote RR Lyrae 
stars, and the photometric errors are shown. 
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Fig. 2. — The differences in fluxes between two synthetic spectra of RGB stars at the HB level, (a, 
b) For the observed variations in lighter elements (A[N/Fe] « 0.47, A[Na/Fe] « 0.53, A[A1/Fe] ss 
0.20, and A[0/Fe] « —0.45). (c, d) CNO and Na are enhanced by 0.3 dex, while heavier elements 
are fixed, (e, f) Elements heavier than Al are enhanced by 0.3 dex, while lighter elements are fixed. 
Panels (g) & (h) are normalized responses for CTIO-4m U and I filters, while the dashed lines are 
for the Johnson filters (see text). 
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Fig. 3. — Same as Figure 1, but zoomed around the SGB and RGB regions in U vs. U — I CMD. 
Denoted by blue and re d dots are, respectively, "Calcium-normal" and "Calcium-rich" stars from 
J.-W. Lee et al.l (|2009al ). Note that two discrete RGB sequences are well separated in calcium 
abundance. 
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Fig. 4. — Our population models for NGC 1851 on the observed C MP in U vs. U — I. The 
models in insets are for V vs. V — I, which should be compared with iMilone et al.1 (|2008l ) for the 
observed SGB split in 7 vs. V — I CMD. Panel (a) is for the case in which only metallicity is 
different (AZ = 0.0004) between the two subpopulations, while panel (b) is for the case where both 
helium and metallicity are enhanced in Pop-2 as in Table 1. The red lines & dots are for more 
metal-rich Pop-2, while the blue lines & dots are for Pop-1. The dashed line in panel (a) shows the 
line blanketing effect of lighter elements only (see text). The crosses denote RR Lyrae stars. The 
apparent distance modulus of (m — M)f / = 15 .61 is adopted, and the color- T e g- transformation and 



bolometric correction from lGreen et al.1 (|1987l ) was used. 
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Table 1: Input parameters adopted in our best simulation of NGC 1851 



Parameter 


Population 1 


Population 2 


Z 


0.0012 


0.0016 


Y 


0.2324 


0.282 


[a/Fe] 


0.3 


0.3 


Age 


10.7 Gyr 


10.6 Gyr 


rf 


0.59 


0.59 


AM h 


0.2337 


0.2218 


<Jm c 


0.015 


0.015 


Population ratio 


0.7 


0.3 



" IReimera (|1977f ) mass-loss parameter 
6 Mean mass-loss on the RGB (M Q ) 
c Mass dispersion on the HB (M ) 



